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C/3 ' Abstract 
O 

We propose a four-level model where sub-Doppler resolution as well as enhanced absorption of a weak 
probe field are realized by using two coherently driving fields. We show that spectral resolution can be 
Q ' improved by modifying the coherent fields intensity and frequencies. 

^ PACS number(s): 42.50.Hz, 42.50.Ct, 42.50.Gy 
• ^ I. 

I Introduction 

1 I ' Multi-photon transitions are embedded in a variety of optical effects based on correlated absorption and emis- 

sion of two or more photons. Step- wise and multi- photon processes can be distinguished by their frequency- 
' correlation properties Q |^. Sub-Doppler resolution for an inhomogeneously broadened medium based on 
multi-photon processes has attracted much attention in recent decades. The physical essence for this method 
\^ lies in the fact that a detuning for multi-photon transition stipulated by Doppler frequency shift can be reduced 
, or eliminated by adopting appropriate light propagation direction, since the detuning is the sum or difference of 
multiple single-photon transition detunings. Sub-Doppler nonlinear optical resonance as appearance of quantum 
coherence and interference in the context of frequency-correlation properties of the coherent components have 
been proposed in |^ . A widely used method is Doppler- free two- or multi-photon absorption where the velocity 
dependent detuning could be removed in the case that the wave vectors of the interacting beams sum down to 
zero (for a review, see [Q). This type of sub-Doppler spectrum is characterized by a large detuning from the in- 
termediate resonance, measuring the fluorescence from the upper level as well as by small cross sections for these 
multi-photon processes. Enhanced by intermediate resonance sub-Doppler processes based on the use of strong 
driving fields and on a change of frequency-correlation properties of resonant multi-photon process in strong 
Q-i' resonant fields, have been proposed in |^ and further developed in j5[ |^ (also in 1^). Substantial enhancement 
^ , in absorption, gain and fluorescence controlled with the auxiliary appropriately propagating electro-magnetic 
wave has been predicted. An interference nature of the spectrum modification has been stressed, which implies 
that along the growth of the absorption (gain) in certain spectral intervals, integral over the frequency may even 
■ decrease. A role of increase of intensity of a probe field as well as features of Doppler-free lasers were explored 
too. 

Control of atomic response with intense coupling lasers has been a subject of many intensive studies in the 
context of electromagnetically induced transparency (EIT) , amplification without inversion (AWI)0, H, 

enhancement of the refractive index without absorption |Q and so on (for review see ^ |13[|). Sub- 
Dopplcr resolution by using intense coherently driving field(s) at other transition(s) has been recently further 
explored in In Ref. probe weak field absorption spectrum in three-level schemes was considered. It 

was pointed out that the linewidth of one of the two Autler-Townes absorption peaks can be reduced by match 
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of the coupling field intensity and frequency. In a later work 1 15 , an experimental observation of sub-Doppler 



linewidth in a Doppler-broadened A-type Rb atomic system was reported. Related works to reduce Doppler 
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broadening with atomic coherence effects could be found in Ref. [g^ and multiple sub-Doppler lines have been 
shown to achieve with a strong coupling field and a saturation effects in a three-level system [pTj. 

In the recent papers ||l8|, cancellation of Doppler broadening by applying different frequency fields was 
proposed with one or two coherently driving lasers. As it has been outlined, sub-Doppler structures can be 
induced without population redistribution. At that not only absorption and gain, but four-wave mixing output 
can be enhanced through Doppler-free coupling controlled with auxiliary co- or counter-propagating driving 
electromagnetic radiation. 

In this paper, we propose a four-level scheme, where both upper and lower levels of a probe transition are 
coupled to other levels by strong coherent fields. In this scheme, four absorption peaks could be found on 
account of the fact that both the upper and lower states are split into Autler-Townes doublets. Sub-Doppler 
resolution is achievable because the modified two-photon transition occurs contributing in the probe process 
and two-photon detunings, resulting from the Doppler shift, could be reduced by choosing appropriate optical 
geometry. A crucial role of the ratio of the frequencies of the coupled transitions is outlined. The results 
obtained on the basis of the developed theory are accompanied with numerical simulation addressed to real 
experimental schemes. Despite the decreased integral intensity, we predict enhanced laser-induced sub-Doppler 
absorption peaks along with transparency windows. The features are interpreted in the terms of quantum 
coherence and interference processes. 

The paper is organized into four sections. In the next section, we present a four-level model and the density- 
matrix equations describing the system. The solution in the limit of a weak probe and steady-state condition is 
obtained and the absorption spectrum is then analyzed, first without Doppler broadening. In Sec. Ill, taking 
into consideration Doppler broadening, we demonstrate that the multiple sub-Doppler lines as well as enhanced 
absorption at the resonance could be obtained in this scheme. From numerical illustrations of the effect in 
a practical medium, the conditions for and features of Doppler-free resonance, stipulated by compensation of 
Doppler shifts with light shifts are discussed. In Sec. IV, we summarize the results. 
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Model and absorption spectrum at homogeneously broadened 
probe transition 
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Fi gUre 1 '. The energy level schemes 
under consideration. 



We consider a closed four- level scheme shown in Fig.|l|(a). In this scheme 
two coherent driving fields Ec and Ed with coupling Rabi frequencies Qc = 
— Ec • di3/2?i and fid = —Ed • d24/2?i interact with the transitions labeled 
|4) — |2) and |3) — respectively. The transition |4) — |1) is probed by 
the weak nonperturbating field Ep. Absorption index ap for this field, re- 
duced by it's value a^g at ojp — UJ41, flc — fid — found as ap/a^o = 
-^^{Xp/Xpolj where Xp is corresponding dressed susceptibility. The later is 
convenient to calculate with aid of dressed nonlinear polarization and 
density matrix pij as P^^{u)p) = NxpEp/2 = Npijdji {dji is transition 
electrodipole moment, N - number density of atoms). In the framework of 
semi- classical theory and using the standard density matrix formalism with 



the rotating wave approximation, the description equations of this scheme can be written in general form as: 



L. 



where Lij = d/dt 



nnPnn 

r • 



qn - 'i[V,p]nn + '}2m>n^rnnPmm] LliPli = LpPp -i[V, p]u (etc.), 



(1) 



Vii — flp ■ expjiApi}; flp = — Ep • di4/2?i; Ap — ujp — ujh (etc.) are frequency de- 



tunings from the corresponding resonance; Tij - homogeneous half- widths of transitions (in absence of collisions 
Tmn = i^m + ^n) f'^)] = 7"j " invcrse lifetimes of levels; "/mn - rate of relaxation from the level m to n, 



qn 



rate of incoherent exitation to a state n from the underlying levels. 



We represent density matrix elements as pi4 = rp ■ exp{iApt}, P33 = etc. Then in a steady-state regime a 
set of density-matrix equation may be reduced to the set of algebraic equations for the amplitudes . Analytical 
solution of this system of 10 coupled equation both for closed and open four-level scheme is given in |2C[| . 
In the case under consideration incoherent exitation of the upper levels as well as relaxation from the level 3 to 
2 are supposed to be negligible small. Then ^2 = ^4 = Tc = and equation system (|^) reduces to 

PpTp = iflpTi + iri2flc - ifldrsi, ^2^12 = irpfl^ - ifldr32, (2) 
-P347'34 = -ifl^p + ir^ftp + ir32flc, P23r23 = -iflcr43 + ir2ifld- (3) 
PdTd = ifldin - rs), Fara = -2Re{inyd}, n ^ I - rg. (4) 
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-iAp,d, P12 = ri2 + i(Ap-Ac), P34 =r34 + i(Ap-Ad), P23 = r23 + «(Ae- Ap + A^). 



For atom moving with speed v, Doppler shift of resonances must be taken into account by substituting Aj for 



Under considered approximation sohition for two level system |1) 



kjV. 
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^dPd 



T2(l + ac,) + A2' r2(l + ac) + A2' T.Td 

Dressed susceptibility for the probe field is found in the notations of pi], |20|] as: 



1 3) is found apart of the other elements. 

(5) 
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Ti (1 + 9h + Vb) - jri - r3)(l + gs - ^6)91 
1 + 55 + I's + 54(1 +95- ve) + V4{1 + W5 - 56) 
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At fid — all gi ~ and the equation (|^) 
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reduces to that describing A scheme: 
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i{Ap 



(rp + iAp)[ri2 + z(Ap- Ae 
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rp[ri2 + z(Ap- A,)] 
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Ac + ?:(ri2 + rp 



Ac + iiTu-Tp) 
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Following to , we introduce frequency-correlation factor 
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The denominator in displays two resonances. At — > we obtain 5i Q 
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iri2. This 



the HWHM is Fp, which corresponds to onc-photon resonance with no 



indicates one resonance at Ap — 0, 
correlation with ojc {Mi = dSi/dAc = 0). The second resonance at Ap = Ac is of HWHM F12 that corresponds 
to two-photon resonance, fully correlated with ujc {M2 = d52/dAc = 1). With growth of the coupling Rabi 
frequency flc the resonance becomes split in two component, their frequency-correlation properties modify, and 
Ml « M2 ~ 1/2 at \flc\^ ^ {|Acp,rp^2}, which neither correspond to one- nor to two-photon processes. 
HWHM of the resonances becomes also nearly equal to each other and to (ri2 +rp)/2. Note that all the effects 
are determined by the coherence pi2 induced in the transition |1) — |2) by two coupled fields, and that always 
Ml + M2 — 1. More detailed discussion can be found in ||l|, ^, 0, [l^ . 

In the alternative case flc — G all Vi — Q and the equation (^) reduces to that describing V scheme: 



Qp 
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[r34 + i{Ap - Ad)] - (ri - r3)|r>rf| V(Fp + iAp) 
{Vp + iAp)[T:i4 + i{Ap-Ad)] + \fldV 



(10) 



The structure of the denominator is similar to (0) and determined by the coherence p^^. Additional term in 
the nominator is stipulated by the coherence induced at the transition |1) — |3) with not zero populations of 
the levels unlike the transition |2) — |4). Alongside with the coherence P34 this term is a source of the nonlinear 
interference effects (NIEF) in absorption (gain) and refraction. Specific features of NIEF in coupled Doppler 
broadened A, V and ladder schemes were explored in |l^. It has been outlined that frequency integrated 
absorption index is proportional to ri — only, i.e. it's change is determined only by the population change. 
As it was first outlined in ||, 19 1, indeed NIEF give rise to difference in the line shapes of pure absorption 



and emission spectra. As the consequence of this effect, the appearance of amplification without inversion on 
the base of NIEF has been predicted. Amplification without inversion was introduced and its features were 
explicitly analyzed and illustrated for the model of neon transitions in the early publication 

Specific feature of the case under consideration in this paper is that both of the levels of the probe transition 
can be driven independently. This gives rise to multiple resonance structure, corresponding to the roots of 
denominator in (^. If resonance splitting is much greater than their widths, we can set Vij — 0, and the 
resonance values of Ap are described by the equation: 



[Ap -Ad- Ac][Ap(Ap - Ad)(Ap - A^) - (Ap - Ac)fll - (Ap - Ad)fll 
-Ap[(Ap - Ad)fll + (Ap Ae)f^2] + {nl - fll)^ = 0. 



(11) 
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Basically, this equation possesses four roots, that indicates appearance of four nonlinear resonances, determined 
by splitting of each of the levels |1) and |4) into two quasi levels. 

In the case = —Ac = A resonance positions are given by the equation: 

Al = (A2 + 2nl + 2n'J/2 ± V[(A2 + 2f]^ + 2J72)/2]2-(f)^-f]2)2. (12) 

At = \^c\ = \^\ two resonances merge in one not shifted resonance at Ap^^ — 0, the other two resonances 
are given by 

A^''^^ = ±VA2+4r!2. (13) 
In the case A^ = Ac = resonance positions arc found as: 

A^i'') = ±{\nd\ - l^icl); a'^'^^) = ±{\nd\ + m). (14) 




Figure 2: Absorption spectrum at homogeneously broadened probe transition |1) — |4) controlled with two driving field. The 
detunings and Rabi frequencies of the driving fields are: Ac = A^j = 0. (a) |f7c| = SFp. (b) = SFp. (c) Ap = 0. 




Figure 3: Absorption index at homogeneously broadened 
probe transition |1) — |4) controlled with driving fields. 
Ad = Ac = 0; 1 - (dash-dot) = 120 MHz, Uc = 200 
MHz; 2 - (solid) = Qc = 120 MHz; 3 - (dash) = 0, 
Qc = 120 MHz; 4 - (dot) = 120 MHz, He = 0; 5 - 
(solid) = Cc = 0. 



We shall illustrate major outcomes of the paper with 
numerical simulations addressed to the conditions of the 
experiment Transitions |1) - |3) - |2) - |4) - |1) in 

Fig. |l|a are attributed to those of the sodium dimers Na2: 
X'T,+ {v" = 0, J" = 45) - A'E+(6,45)(Ad = 655 nm) - 
Xil]+(14,45)(A32 = 756nm)-Bin„(5,45)(Ac = 532 nm)- 
X'I]+(0,45)(Ap = 480 nm). The Doppler widths of the 
transition at wavelength Ap = 480 nm at the temperature 
about 450 C is 2Dp « 1.7 GHz. Then the Boltzmann's 
population of the level |2) makes about 1.5% from that of 
the level |1). The following relaxation parameters are used: 
r4 = Fg = 120, Fs = Fi = 20, 742 = 5, 741 = 10, 732 = 4, 
731 = 7, F12 = 20, F34 = 120, F23 = Fc - Frf - Fp = 70, 
(all in 10^ s~^). For numerical simulation we have used the 



full set of the equations from ||T^, ^ , accounting for various 
relaxation transitions and not zero population of the level 
|2). 

We begin our discussion with the case when the Doppler-broadening is not included. Performing numerical 
calculation for Eq. (6), we depict the absorption profiles in Fig. ||. From Fig. || a and b, it is easy to find the 
typical EIT and absorption spectrum with {fldl = (Fig. ||a) and \flc\ = (Fig. |^). Under increasing Rabi 
frequency |llc| (|f^d|), the four-peak spectrum forms, when the driving fields become comparable to or larger 
than the spectral width Fp of the atomic transition. The emergence of the four-peak spectrum is attributed to 
the fact that levels |1) and |4) have been driven into two dressed states (Fig. p. As {fldl approaches \i^c\, the 
four-peak spectrum degenerates into a three-peak spectrum since two dressed-state transitions have the same 
resonant frequency and contribute to the same central component. As already known, a three-peak spectrum 
occurs in resonance fluorescence (RF) in the strong field limit of a two-level system (for a review, see ||]). Even 
though the dressed-state diagrams are similar in the two cases, the four-peak spectrum can not appear in RF 
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because in the case of RF the dynamic Stark sphtting for the two doublets must be the same, while in our system 
the two doublet splittings are controlled separately by different driving field intensity and frequency. We would 
like to stress that the multiple EIT windows can appear simultaneously with the multiple absorption peaks. In 
Fig. I^c, the absorption amplitude A at Ap = is plotted against the Rabi frequencies for two resonant driving 
fields. An important feature is that the maximum absorption occurs under the condition that the two Rabi 
frequencies are equal, which is due to the fact that the two dressed-state transitions simultaneously contribute 
to the central component. This result may provide a way to measure the field intensity as well as atomic 
parameters. 

Fig.^ depicts absorption spectrum for the molecules with the velocity projection on the propagation direction 
of the coupled fields v = 0. The plots illustrate dependencies, described by the equations (|l^) - The 
appearance of the four-peak spectrum is attributed to the fact that levels |1) and |4) have been driven into two 
dressed states. Direct computing shows that frequency integral absorption for the plot 3 is the same as for 5, 
whereas for the other plots it is about 0.5 of that for the graph 5. This is due to strong saturation of population 
difference at the transition |1) - |3). 

Ill Coherence induced resonances in Doppler-broadened medium of 
sodium dimers 

In the previous sections, Doppler broadening has not been accounted for. In order to consider this effect, 
Doppler shifts must be introduced in the three interacting field resonance detunings in Eq. Below we will 
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Figure 4: Absorption spectrum at inhomogeneously broad- 
ened probe transition |1) — |4) (the scheme Fig|l](a)) con- 
trolled with two driving field. All fields are co-propagating. 
Ajj = Ac = 0, He = 1 GHz, is as follows: 1 — (solid) = 1 
GHz; 2 - (dash) = 0.7 GHz; 3 - (dot) = 0.3 GHz; 4 - 
(solid)- Cc = = 0. 



Figure 5: Compensation of residual Doppler broadening and 
sub-Doppler structures at the probe transition |1) — |4) (the 
scheme Fig.h|b) coherently driven with two strong fields. Main 
plot — all fields are co-propagating, 1,-3 — all parameters 
and notations are the same as in Fig.W. 4 — (dash-dot): fl^ = 
0.075 GHz. Insets (a) and (b) are corresponding zoomed peaks. 
Inset (c) — same as plot 1 but: 1 — only field c is counter 
propagating; 2 — only field d is counter propagating; 3 — both 
driving fields are counter-propagating in respect to the probe 
field. 



proceed with analysis of the absorption spectrum under effects of Doppler-broadening. Figure ^ displays same 
dependencies, but for velocity averaged absorption index and co-propagating probe and driving fields. While 
number of peaks and their positions are identical to the previous case, the specific feature of the inhomogeneous 
broadening is that plot 1 displays enhanced absorption in the central component despite the fact that integral 
intensity of this graph is about 2 times less than that for the plot 4 (due to the strong saturation of the 
populations at the transition |1) - |3)). In absence of driving fields only small fraction of the molecules (about 
Tp/Dp « 10~^) can be coupled by probe field. Overlap of the two dressed transitions gives rise to increased 
amount of the molecules coupled with the probe field and consequently to enhanced absorption as compared to 
not perturbed absorption in the center of Doppler broadened transition. 
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As it was shown in resonance requirements for molecules at velocity v and at maa;{|ridp, A^P} ^ kdU 
{u is thermal velocity, kd - wave number) are described by the equation 

Ap - kpV = 61^2 - Mi,2kdV. (15) 

The equation indicates that cancellation of Doppler shifts is possible even at kp ^ kd, if kp < kd- This is due to 
the interplay of the Doppler and ac Stark shifts giving rise to variation of Mi in the interval 0...1/2, whereas M2 
- in the interval 1...1/2. Possible enhancements of cross-section of optical processes through concurrent coupling 
of molecules from wide velocity interval by means of the above outlined approach are shown in Fig. |^. The 
frequencies of the probe and Ed fields are interchanged (see Fig. |l|b), all other relaxation parameters remain 
the same. Plot 1 displays enhanced sub-Doppler resonances with the FWHM comparable with the natural 
linewidth (see insets (a),(b))for co-propagating wave. Dramatic change as compared to similar graphs in Fig. ^ 
is explicitly seen. The other plots show evolution of the structures with variation of the intensity of one of 
the driving field. As plot 4 indicates, the sub-Doppler structures exist only under certain ratio of the driving 
field intensities. Even more dramatic change occurs while driving fields (especially shorter wavelength one) are 
counter propagating (see inset (c)). 

IV Conclusion 

In conclusion, we have explored spectral properties of the molecular transitions in the case that both upper 
and lower levels are coupled to other atomic levels and Doppler effects play an important role. The features 
of the four-peak or three-peak spectra induced by the driving fields are investigated. Despite the dominated 
Doppler-broadening, two or more peaks may possess the sub-Doppler resolusion. A crucial role of the ratio 
of the frequencies of the coupled transition is shown. With manipulation the detunings and/or intensities of 
two coupling fields, one can improve the spectral resolution. The enhanced absorption in the sub-Doppler peak 
is shown to be realized, whereas the peak can be induced in the center of the inhomogeneously broadened 
transition. The predicted effects are attributed to the quantum coherence and interference processes, while 
frequency-correlation properties of multi-photon processes experience substantial modification with the growth 
of the driving fields intensities, which leads to corresponding dramatic changes of the role of Doppler effects. 

The similar sub-Doppler technique can also be implemented for other atomic coherence effects, such as EIT, 
LWI and for enhancing various resonant nonlinear optical phenomena in inhomogeneously broadened media. In 
one of our previous papers [ p2| , the inversionless amplification has been considered with the important feature 
that incoherent excitation to the upper level is not necessary in the same scheme as considered in this letter. 
Similar scheme is often used in quantum optics. Most recently, electromagnetically induced absorption was 
studied in a similar four-level system We believe that, by utilizing this sub-Doppler technique, the above 
effects in inhomengeneously broadened media can be further enhanced and manipulated. 

We want to stress that in optically thick media the processes similar to optical parametric amplification and 
involved in the schemes Fi g. |l| a,b may play a crucial role, that imposes dramatic consequences on the features 
of the output probe signal pif . 
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